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Abstract 

We show that the quantum properties of some Josephson SQUID devices are 
described by a boundary sine Gordon model. Our approach naturally describes 
multi-junction SQUID devices and, when applied to a single junction SQUID (the 
rf-SQUID), it reproduces the known results of Glazman and Hekking. We provide a 
detailed analysis of the regimes accessible to an rf-SQUID and to a two-Josephson 
junction SQUID device (the dc-SQUID). We then compute the normal component 
of the current-response of a SQUID device to an externally applied voltage and show 
that the equation describing the current-voltage characteristic function reduces to 
well-known results when the infrared cutoff is suitably chosen. Our approach helps 
in establishing new and interesting connections between superconducting devices, 
quantum brownian motion, fermionic quantum wires and, more generally, quantum 
impurity problems. 
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1 Introduction 



Quantum effects in Josephson junctions have been by now investigated for 
quite a long time, both experimentally [T] and theoretically [2J: the most 
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important manifestation of quantum fluctuations in these systems being the 
macroscopic quantum tunneling of the phase across a current biased junction 
and the consequent observation of quantum phase slips events at a bias rather 
close to the Josephson critical current pQ. 

Early studies of quantum phase fluctuations in Josephson systems have been 
carried using either the Caldeira- Leggett [3] or the electromagnetic environ- 
ment [1] models: within this framework, it has been shown that macroscopic 
quantum tunneling causes a finite voltage to appear for any finite current, lead- 
ing to a non-linear I-V characteristic which, in the limit of zero temperature, 
reads as V = AP with 7 depending on the impedance of the superconducting 
leads and A being proportional to the bare tunnel matrix element between 
two adjacent minima of the potential [5]. The dual result for a voltage biased 
junction shows that the dc-current in such a junction is proportional to the 
square of the unrenormalized Josephson energy [BJ. 

As already pointed out in Ref.[7], in the above mentioned approaches the ef- 
fective boundary conditions for the quantum fluctuations of the "environment 
modes" do not depend on the Josephson energy of the junction. While this as- 
sumption is perfectly legitimate for weak fluctuations of the phase of the order 
parameter across the junction, better care should be used if these fluctuations 
are strong as it may well happen for one- dimensional superconductors, where 
phase fluctuations diverge logarithmically with the length of the system, open- 
ing the undesirable option that the Josephson energy could average to zero 
for diverging random fluctuations of the phase. 

In their seminal paper Glazman and Hekking [7J showed that a finite renor- 
malized Josephson energy may arise since the junction itself affects the fluc- 
tuations of the environment. They considered an rf-SQUID, i.e. a thin super- 
conducting loop containing a Josephson junction, phase-biased by threading 
a flux $ ext through the loop. The environment modes of the loop consist of 
plasmon modes propagating through a superconducting wire with a soundlike 
dispersion law [8]. Their analysis was based on a perturbative renormalization 
group approach exploiting the similarity of the effective action, emerging from 
integrating out the fluctuations "away" from the junction, with the one used 
to analyze the quantum brownian motion in a periodic potential [9]. They 
found how the Josephson coupling Ej is renormalized when high energy de- 
grees of freedom are integrated out and showed that the renormalization leads 
to values of the Josephson critical current which are smaller than one would 
expect from the mean field theory result J° = 7rAG/2e, where G is the con- 
ductance of the junction and A is the superconducting gap in the loop. The 
RG approach is supplemented by an instanton analysis [10] needed to compute 
the effect of the macroscopic quantum tunneling on the phase dependence of 
the Josephson current for relatively large loops. 
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Generalized SQUID devices are much studied today mainly because of their 
relevance for the implementation of flux qu-bits [11]: it is by now widely ac- 
cepted that multi-junction SQUID devices [12] are the most promising real- 
ization of a flux qu-bit due to their reduced size, which renders them less 
sensitive to decoherence effects [13] induced by the interaction with the envi- 
ronment. As it is well known, an rf-SQUIDs [H] may be represented as a two 
level system only when the superconducting loop is rather large and, thus, 
unfortunately more sensitive to interactions with the environment. 

In this paper we evidence that all the non-perturbative renormalization effects 
induced by quantum fluctuations of the phase of the superconducting order pa- 
rameter [7J are described by the two-boundary sine Gordon model [H] , which 
becomes then the pertinent description of the dynamics of the fluctuation 
field and of a Josephson junction with arbitrary nominal strength Ej. As we 
shall see, the quantum behaviors of rf-SQUIDs and of multi-junction SQUID 
devices, may be analyzed within the framework of boundary field theories |16j . 

A two boundary sine Gordon model description of SQUID devices is possible 
since the dynamics of the fluctuations of the phase of the superconducting 
wire may be described as a spinless Tomonaga-Luttinger liquid (TLL) [TTJ . In 
fact, for a superconducting wire of length 2L and cross-sectional area S, the 
Lagrangian describing the phase fluctuations of the order parameter may be 
written as [TJ: 




where n s is the superfluid density, m is the electron mass, l/e c is the charac- 
teristic inverse charging energy per unit length of the loop: 

I = ! (2) 



where e is the dielectric constant of the medium surrounding the loop, while R 
is the distance from a metallic screen, which, together with a, fixes the lead's 
capacitance. 



Upon setting h = 1, if one defines the "Luttinger" parameter g = ^\Jn s S/ (me c ) 

and the plasmon velocity v p \ = ^Jn s Se c / (4m), one immediately sees that the 
Lagrangian (CQ) may be written as the one describing a spinless Luttinger liquid 
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1181: 




(3) 



2 

Physically, 1/g is the dimensionless zero-frequency impedance of the super- 
conducting wire [TPU] . Due to its dependence on S, g may take values in a 
rather wide interval of real numbers: for g > 1, Eq.([3]) describes a TLL with 
an attractive interaction, while, for g < 1 the interaction is repulsive. The 
value g = 1 describes the very special case of an essentially free theory. Differ- 
ent superconducting wires may realize different values of g: for instance, for 
rather clean aluminum wires, one can easily attain values of g such that g > 1 
and, thus, realize a wire for which the dynamics of the phase fluctuations is 
described by an attractive TLL. [TJ. 

For g > 1, the two boundary sine Gordon model description of a SQUID device 
allows to make very general statements regarding the regimes accessible to the 
system for a finite loop-size L. Namely, there will be a perturbative (in Ej) 
weak coupling regime accessible when the relevant Josephson couplings are 
small and a non-perturbative strong coupling regime accessible when Ej be- 
comes large; most importantly, there will be a renormalization group invariant 
length scale such that for L < the SQUID device is in the perturba- 
tive weak coupling regime, while it is in the non-perturbative strong coupling 
regime for L > L*. Intuitively speaking, is (2tt times) the radius for which 
the ratio between the Josephson energy and the magnetic energy due to the 
loop self-inductance equals one; the existence of an "healing length" is 
generic to a boundary sine Gordon model describing a TLL [21]. At variance, 
for g < 1, the theory lies within the perturbative regime at any L. Therefore, 
there is no healing length in this case and no universal behavior is then at- 
tainable. In the following, we shall show that these properties are shared by a 
variety of SQUID devices. 

We investigate also the stability of the two (weakly and strongly coupled) 
regimes by looking at the relevance/irrelevance of the operators describing 
the boundary interactions. It turns out that these operators are either vertex 
operators or instantons (phase slips) of the plasmon field $, depending on if 
the system is within the weakly, or the strongly coupled regime. As evidenced 
in the following, vertex (instantons) are a relevant (irrelevant) perturbation 



2 Our normalization of g renders the boundary interaction marginal when g = 1. 
It differs by the one used in Ref.[T9] by a factor of 1/2. 



4 



of the weakly (strongly) coupled regime for g > 1, while they are irrelevant 
(relevant) for g < 1. 

It emerges the picture that, for g > 1, it always exists an healing length 
L* such that, for L > L*, the Josephson coupling affects the dynamics of 
the plasmon field as a boundary condition: in this regime it is impossible to 
disentagle the Josephson junction from the environmental modes. For g < 1, 
instead, this situation is never attainable and the Josephson coupling has to 
be always regarded as a quantum impurity embedded in a superconducting 
wire. 

The strongly and weakly coupled regimes may be probed by looking either at 
the dc-conductance, or at the behavior of the Josephson current as a function 
of the magnetic flux treading the superconducting loop. The two-boundary 
sine Gordon approach provides a systematic and powerful method to evaluate 
the relevant currents. In fact, the correlators of the plasmon fluctuation $ 
field(s) may be exactly computed and this allows for a systematic derivation 
of the leading corrections to the Josephson current for all values of the applied 
flux ip = $ cxt /$o ( $g = h/(2e)), which, for cp = ir, yields the results of the 
instanton analysis carried in Ref. [7|; in addition, we shall show that it also 
allows for a systematic derivation of the response of a SQUID device to an 
applied external voltage. Since applying a finite voltage V introduces a new 
energy scale e*V, with e* = 2e, the dc-current response of the system strongly 
depends on the ratio between the new energy scale, and the intrinsic scale 
provided by the level spacing s = 2nv p \/L of the plasmon field(s) $. Our 
computation of the dc-current clearly shows that, when e*V > s, one gets 
the well known Kane-Fisher formula [9], while, for s > e*V, one gets the 
results obtained by Glazman and Hekking [7]; in addition, when e*V > s, the 
dc-current dominates over the Josephson current. 

Boundary field theories have by now become relevant in several different con- 
texts. In condensed matter theory , they are mostly generalizations of quantum 
impurity models, which may be described by using the TLL paradigm [T7]; for 
instance, boundary interactions appear in the analysis of the Kondo problem 
[22], in the study of a one-dimensional conductor in presence of an impurity 
[23], in the derivation of tunneling between edge states of a Hall bar [21], in 
the study of quantum brownian motion in a periodic potential [25] and in 
the analysis of the phases accessible to networks of quantum wires [19J and 
Josephson junctions [26] , [21] . The TLL paradigm evidences that many inter- 
actions are simply diagonalizable in the basis of pertinent collective bosonic 
modes, and that non-diagonalizable interactions correspond to exactly solv- 
able Hamiltonians, such as sine Gordon models [2Z] • Recently, boundary field 
theories have been investigated in the context of string theories. For instance, 
in studying tachyon instabilities [28] , one is faced with the fact that the space 
of interacting string theory [22] is mapped onto the space of boundary per- 
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turbations of conformal field theories [30], and that the renormalization group 
flow determined by boundary perturbations may be identified with tachyon 
condensation [3T] . Affleck and Ludwig [H2] showed that the boundary entropy 
g is decreasing along the renormalization group trajectories, triggered by the 
boundary interactions. It is remarkable that this entropy has been proposed 
as a measure of entanglement for boundary states |32j . 

The paper is organized as follows: 

Section 2 shows that quantum behaviors of an rf-SQUID are described by 
a two-boundary sine Gordon model. There, we shall analyze the weak and 
the strong coupling regime accessible to the device for a finite size of the 
superconducting loop. We show that, for g > 1, there is a renormalization 
group invariant scale - the "healing length" L*- [7\ marking the crossover 
between the weakly coupled and the strongly coupled regime. 

We first analyze the weakly coupled fixed point corresponding to mixed bound- 
ary conditions (i.e. Dirichlet at the inner boundary and Neumann at the outer 
boundary) and we determine the scale dimension of the pertinent boundary 
operators and the scaling equations for the running coupling constants using 
renormalization group (RG) methods; this allows for a perturbative computa- 
tion of the beta-function for all the couplings. Computation of the canonical 
ground state energy shows that- at the weakly coupled fixed point- the Joseph- 
son current has a sine dependence on the flux treading the loop. Then, we 
analyze the strong coupling fixed point, corresponding to Dirichlet boundary 
conditions at the outer and inner boundaries of the rf-SQUID. Evaluation of 
the canonical ground state energy shows that- at the strongly coupled Dirich- 
let fixed point- the Josephson current exhibits a saw-tooth dependence on the 
external flux. 

Due to the zero modes, appearing in the strong coupling regime as a conse- 
quence of the Dirichlet boundary conditions, the ground state is doubly de- 
generate when (p = 7T + 2ixk: when these levels are confined away from the rest 
of the spectrum, the rf-SQUID realizes a qu-bit (See, for instance, Ref.JTT]). 
Our analysis points out that this realization is possible only for g > 1 and for 
sizes of the superconducting loop bigger than the healing length, and, thus, 
for a rather big self-inductance of the device. Finally, we analyze the effects 
of quantum phase slips on the ground state degeneracy of the rf-SQUID at 
the strong coupling fixed point: an explicit instanton computation shows that, 
for g > 1, the degeneracy is lifted by an amount proportional to the instan- 
ton fugacity, and that the instantons smoothen the edges of the saw-tooth 
phase-current relation. 

In section 3 we analyze a SQUID-device involving only two Josephson junc- 
tions, of arbitrary strength Ej t %, E Jj2 (the dc-SQUID). Under the assumption 
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that the two plasmon fields living in the two branches of the superconducting 
loop (see fig.l) are described by TLLs corresponding to the same g, we show 
that also the dc-SQUID is described by a two-fields generalization of the two 
boundary sine Gordon model. We show that also here, for g > 1, there is an 
unstable weak coupling fixed point and a stable one at strong coupling and 
that the crossover from one to the other is characterized by an healing length 
L*, whose value depends, of course, only on the smallest of the two Josephson 
couplings. When g < 1, the only stable fixed point is at weak coupling. 

In section 4 we analyze the normal component of the current-response of the 
rf- and dc-SQUIDs to an applied external voltage. There, we evidence that 
a computation of the time-independent current using the $-field(s) correla- 
tors derived from the two-boundary sine Gordon model yields the well-known 
Kane-Fisher formula when the infrared cutoff e*V > s |9J, as well as the 
Glazman-Hekking result [7], when the infrared cutoff if e*V < s. It is remark- 
able that, as a result of the boundary field theory approach to superconducting 
SQUID devices developed in this paper, one gets a current-voltage equation 
valid for any choice of the infrared cutoff. Furthermore, for a dc-SQUID, the 
boundary field theory approach evidences the possibility of new remarkable 
interference effects (even in the absence of an external applied magnetic flux 
(!)) between the currents flowing across the two junctions. 

Section 5 is devoted to our conclusions and final remarks. 

In order to be pedagogical and self-contained, a few pertinent mathematical 
details are summarized in the appendices. 



2 The rf-SQUID 



From Eq.Q one may easily obtain the Hamiltonian H r f for an rf-SQUID 
(i.e. for a superconducting wire cut by a Josephson junction with nominal 
Josephson coupling Ej and pierced by a magnetic flux ip). By neglecting the 
charging energy of the junction [7J, H r j is given by: 

L 

Hrf = Ho + H B = JL f dx 
Zrr J 

-L 

-£ jC os[$(L)-$(-L)] , (4) 

where H and Hb describe the bulk and boundary Hamiltonian, respectively. 

The flux ip may be accounted for by a redefinition of the field $ as $ — ► 
$ — <pfr- By introducing the symmetric and the antisymmetric phase fields 



+ v p i 



dx 
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$ e (x) = ^=[$(x) + $(-x)] and $ (x) = ^=[$(2;) - with < x < L, 

one gets 

-SjCOs[V2$ (i)+y] -(5) 



The field <& e (x) fully decouples from the interaction term in Eq.(jSJ). The Hamil- 
tonian for the field $ (x) (which will be referred to as $(x), from now on), on 
the other hand, is the Hamiltonian for a two-boundary sine Gordon model |15j . 
with Dirichlet boundary condition at the inner boundary, i.e., $(0, t) = 0, and 
dynamical boundary conditions at the outer boundary. The latter are obtained 
by requiring that the energy functional is conserved, yielding [21] 

9 ^^^ + V2E J smiV2<S>(L,t)+< f ] = . (6) 



For Ej — > 0, one has Neumann boundary conditions at the inner boundary, 
i.e., 9 ^Qx'^ = 0> while, for Ej — > oo, one gets Dirichlet boundary conditions, 
\/2<I>(L, t) = —tp + 2nk, k E Z. As it will be clearer in the following, the effec- 
tive two-boundary sine Gordon model provides the relevant renormalizations 
of an rf-SQUID: namely, it gives the renormalized value (Ej) of the Josephson 
energy and yields the correct functional dependence of the Josephson current, 
as a function of the flux if. 

As we shall see in the following sections, as L — ► oo, there are two relevant fixed 
points, namely, the weakly coupled Neumann fixed point (Ej = 0), and the 
strongly coupled Dirichlet fixed point (Ej = oo). Our analysis will evidence 
the existence of an healing length L* = gnVpi/Ej, providing the size of a loop 
for which the ratio between the Josephson energy and the magnetic energy 
due to the loop self- inductance equals 1, such that, for g > 1, at L — L* the 
device crosses over from the weakly to the strongly coupled regime. 

2. 1 The weakly coupled fixed point 

At the weakly coupled fixed point (Ej — > 0), the field $(x,t) obeys Dirichlet 
boundary conditions at the inner boundary and Neumann boundary condi- 
tions at the outer boundary. As a consequence, <&(x,t) takes the mode expan- 
sion 

Hx,t) = -^Y,—Mk n x)e-^ ut , (7) 
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with [a(n), a(n')} = nS n+n ' !0 , and k n = j [n + |j. 

The ground state of the bosonic theory, |GS), is defined in terms of the oscil- 
latory modes as 

a(n)|GS) = ; V>0 . (8) 



Normal ordering the Josephson energy with respect to |GS), leads to 



H, 



Ej /2na 
~2~ \T~ 



: e 



■itp 



: + :e 



v . 



(9) 



where the column :: denotes, as usual, normal ordering. 

The thermodynamic limit (i.e., large L) is attained once the scaling equations 
for the pertinent running couplings have been obtained using the renormal- 
ization group (RG) approach (see appendix A). 



The RG equations are determined by the requirement that the partition func- 
tion Z is independent of the cutoff. In order to compute Z, one should Wick- 
rotate the field $ and resort to the imaginary-time Feynman path-integral 
formalism. The partition function is, then, given by 



J n$e- So T T [e- SB ] 



(10) 



where T r is the imaginary time ordering operator, 5o is the free Euclidean 
action 
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P 



So= 2^J dT J dX 
o o 



1 /<9$\ td& 
v v \ \9t) pl \dx 



'ID 



and Sb is the boundary interaction action at imaginary times 



S 



E, 



B 



dr 



. e i[V2$(ir)+<p] . _|_ . e -i[V2$(ir)+<p] 



(12) 



with $(ir) = $(L,ir). 

In Eq.f ll2p . the renormalized coupling constant Ej = Ej(2-Ka/L)~9 has been 
introduced. Of course, after this redefinition of the Josephson coupling, there 
is not anymore a cutoff dependence of the interaction operator |33j. 
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The boundary action in Eg. (1121) contains only the vertex operators V±\(t) 



±itp 



exp 



±i\/2$(ir) 



:. Due to the mixed boundary conditions on $, the 
operator <3>(ir) contains no zero modes and, thus, there are no selection rules 
on the Kac-Moody charge carried by the vertex operators. As a consequence, 
the Operator Product Expansions (O.P.E.'s) may generate additional vertices, 
with higher periodicity in ip, given by 



VJt) 



e ±niv : exp 



±iV2n<5>( 



IT 



(13) 



The scaling dimension of V n , h n , is derived from the correlation function 
(T T [V n (r)V^ n (r)}) which, for (3 — > oo, is approximatively given by (see Ap- 
pendix B) 



{T T [V n (r)V. n (r')]) 



tanh 



4L 



T — T 



(14) 



At short imaginary time distances (f p i|r -r'|/L< 1) , Eq. (fl4l becomes 



(r r [K(r)y_„(r')]) 



T — T 



. 2nf 

s 



(15) 



leading to h n = n 2 /g. 

As we shall see in the following, despite the fact that h n depends on n, due to 
nonlinearities, the higher harmonics in Eq. (fl3"|) are irrelevant operators, when 
g < 1, while they become relevant for g > 1. 

The O.P.E. yields the pertinent fusion rules for the vertex operators 
V n (r)V n/ (r') «v 



r 



V n+nl {r') 



(16) 



Introducing the pertinent counterterms, one may easily show that the effective 
boundary action looks as 

S B ,Eff = ~ I dr f^\E n e in ^V n (r) + E n e- m *V„ n (r)} , (17) 
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where the renormalization condition is that all the coupling strengths E n = 
(2ixa/L) 9 E n , with n ^ ±1, vanish at a given reference length scale L Q . 

Setting the running couplings g n as 

9^{^]^°^ , (18) 



from the fusion rules given in Eq.( fl6l) . one is able to derive, in principle, the 
/3-function for all the couplings g n . 

For the sake of simplicity, we shall restrict our attention to the first two rel- 
evant terms, proportional to the running couplings gi,g%- The RG equations 
are, then 

j: 7TTT\ =A(gi»g2) = fl ~ ~ ) 91 + 9x92 , (19) 
d\n(L/L ) \ gj 



and 



dg 2 



d\n(a/ao) 



Mi, 92) 



92 + (91 



(20) 



the flux ip fully decouples from Eqs. (119|20j) and, thus, the scaling near the 
Neumann fixed point is independent of ip. 

For g < 1, gi and gi are both irrelevant. Thus, the Neumann fixed point is sta- 
ble and the theory is perturbative in the boundary interaction with effective, 
size-dependent couplings given by 

91 {L) « g x {L ) (j-j ; g 2 (L) « . (21) 



On the other hand, when g > 1, g\ becomes relevant since it grows with L as 

(i_l) 

9i{L) « gi{Lo) 9 ■ (22) 
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By inserting Eq. (!22|) into the nonlinear term of Eq. (l20l) . one gets 



92{L) 



l + 2/g 



L 

To 



(23) 



where the integration constants have been chosen so as to be consistent with 
the renormalization condition ^2(^0) — 0. 



From Eqs. (l2~3"|) . one sees that the nonlinear term in Eq. ffT^]) makes #2 (L) in- 
crease, as soon as one starts the scaling flow, even when 1 — 4/g < 0. This is 
a remarkable result, since it implies that, as soon as the first harmonics in (p 
sets in, all the others follow, due to nonlinearities in the RG equations. 

At the Neumann fixed point, the Josephson current is perturbative in Ej\ of 
course, as L increases, our results are reliable only for g < 1. As a function of 
the flux (p, I[<f), is given by 



lim — S- In 



c$o /3^°o (3 dip 



J Eff[ 



(24) 



where Z E ff[<p] is the (^-dependent partition function, given by 

ZEffW] = J D$e- So T T [e- Ss ^] . (25) 



To the second order in g n , one gets 



e*E, 



2na 



sin(<^) + ^77T sin ( 2< ^) 
9i{L) 



(26) 



which explicitly gives the contribution of the second harmonics to the Joseph- 
son current. 



For g > 1, Eq. (1231) shows that, for L — > 00, the ratio between the running 
coupling constants of the first two harmonics is given by 



92{L) ^ 
9i(L) \L f 



(27) 



showing the relevance of the boundary interaction, since higher harmonics be- 
come more and more relevant, as L —>■ 00. The scaling equations ( fT8|19|20ir2Tl) 
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cease to be valid at the renormalized healing length for which gi(L*) ~ 1, 
i.e., when 

L = L m = 2Ka(jjf\* : , (28) 



yielding that the renormalized Josephson coupling Ej = Ej(L^) is given by 

(29) 



which is the energy where the level spacing of the plasmon field becomes of 
the order of the Josephson energy, i.e., 2nv p \/L* ~ (2na/L Jf )~iiEj [7|. 



2.2 The strongly coupled fixed point 



In this subsection, we analyze the properties of the rf-SQUID in the limit in 
which Ej/(2nv p i/L) is ^> 1: Hj is now the dominating potential term, and 
the field $(L, i) takes values corresponding to a minimum of the boundary 
energy. Therefore, at both boundaries <&(x) must obey to Dirichlet boundary 
conditions, given by 

$(0,t) = ; y/2$(L, t) = -if + 2irk ; (k G Z) . (30) 



Both boundary conditions are satisfied if one chooses for <&(x,t) the mode 
expansion 

if>(x,t)=-—P-fiY^fa(k n x)e ik » v # t , (31) 
L \9^ n 



where P is the zero mode operator and k n 
the set of possible eigenvalues of P: 



-^ + k 
2vr 



(|n). From Eq.(j30|), one obtains 



(32) 



As a consequence, there will be an infinite set of ground states of the oscilla- 
tory modes, corresponding to the possible eigenvalues of P, and denoted by 

\Pk,{0})- 
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Once H r f has been evaluated for the particular solution given in Eq. (l3TI) . one 
obtains the effective Hamiltonian at the Dirichlet fixed point, Hp, 



H 



D 



TTVpi 



2gP 2 + 2j2®(~ n M n ) 



n=l 



(33) 



from which one gets that the partition function, at the strongly coupled fixed 
point 



Z D [ V } = Tr[e-^ 



1 



(34) 



with 77(2;) = n^°=i(l - x n ), and q = exp \-/3^ 



leading to the well-known sawtooth-like dependence on ip of the Josephson 
current as (3 — > 00 [7j, since 



l dmZ D [<p] 
I[ip\ = - km — cx ip 

P->oc (3 dip 



(35) 



of course, [ip] is the integer part of ip (in units of 2%). 

At strong coupling, the degeneracy among the minima of the boundary term 
is removed by the magnetic energy, E% = (— ^ + k) , selecting only 
one eigenvalue of the zero mode of the field $(x, t), except when ip = it + 2irk, 
where = E^_ v When these two levels are confined away from the others, 
the rf-SQUID may operate as a qubit [11] . For L > L*, however, the oscillatory 
modes of the plasmon field determine a renormalization of the parameters of 
the effective two-level system, as it will be discussed in the following. 



2.3 Renormalization of the Josephson current for the rf-SQUIDs 



At strong coupling and for ip = tt, the ground state is twofold degenerate, 
that is, Ej 0) = Ef>. This de generacy is removed by phase slips (instantons), 
of amplitude V2ti, connecting the ground states. 

In order to describe phase slips within the framework of the two-boundary sine 
Gordon model, one has to introduce the dual field of <&(x, t), Q(x,t), defined 
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by 



d$(x,t) _ v pl dB(x,t) 2gd$(x,t) _ dQ(x,t) 



dx 2g dt v p \ dt 



dx 



(36) 



Due to the duality between $ and 6, the effective Hamiltonian given in Eq.(l33 
may be usefully rewritten as 



H 



1 



D 



An(2g) 



dx 



i foey fde s 



(37) 



From the algebra of the bosonic fields, one easily derives [19] the commutation 
relations 



d$(x,t) 
dx 



e( y ,t) 



dQ{x,t) 
dx 



2m8(x — y) 



(3* 



In order to construct the leading correction to the Dirichlet fixed point, one 
may use the Delayed Evaluation of Boundary Conditions (DEBC) approach, 
introduced in Ref.[19]: for this purpose, one has first to consider the most 
general primary field of the bulk theory, namely, a vertex operator of the form 
V n> n(%,t) =: exp i^<&(x, t) + i^Q(x,t) :, and then to obtain the boundary 

perturbation, (t), by evaluating it at x — L, using the pertinent bound- 
ary conditions for the fields <&(x,t) and Q(x,t). Since <&(x,t) obeys to the 
Dirichlet boundary conditions at x = and at x — L, it will not contribute 
to \Qfi(t). At variance, Q(x,t), obeys to the Neumann boundary conditions 
at both boundaries; namely, one has that 

fle(o, t) = geCM) = Q 

dx dx 



implying the following mode expansion for Q(x,t) 



6(x, t) = 0o + 2 ~^2gP + 2i^g-J2— cos 



Tin 



-x 



(40) 



with 6 = (qn — q£)/y/2 and k n = nn/L. 
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As a consequence of Eq. (|40l) , one finds 

0(ir) = e(L,ir) = + ^2^P + z2^E ^(-1)V-^ (41) 



and thus, the generic boundary perturbation at the strongly coupled fixed 
point, H B , may be written as a linear combination of the vertices Va = V fi 
as 

^B = -E^:e^ 9(L) : . (42) 



At the Dirichlet fixed point, the partition function Zd[<p], including the bound- 
ary interaction in Eq. (|42p . is written as 



Tr 



e-^T r exp<|E / dr K : e l ^ e(lT) : 



(43) 



with 9(zr) = e(L,ir). 

Eq.( |43l) reduces to Eq.( |34l) when all the A n 's are set to zero. 

Following the same steps used in section (2.1), one may derive the scaling di- 
mension of Va, ha, from the large- L- limit of the two point correlation function; 
namely, 

-2gn 2 

(44) 



(T T [F fi (ir)Vl fi (ir / )]) ( o) 



L 



\T — T 



where (. . .)(o) denotes thermal averaging defined by the free Hamiltonian with 
Dirichlet boundary conditions for $(x, t) at both boundaries. 

Eq. (1441) implies that ha = gn 2 and that the leading perturbation at the 
strongly coupled fixed point is given by 

■ e(L) . e(L) 

H B = -Y[[: e~ : + : e~ l ~ :] . (45) 



The physical meaning of the operators may be inferred from the commu- 
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tation relation 



[P,:e 



in@(ir) .1 



1 

2^ 



g^O^r) _ < « e(t) 
ax 



\/2 ' 



(46) 



which shows that Va changes the eigenvalue of P by nj a/2. Thus, the boundary 

fields : e ± ^ 0( - lr ' > : describe instanton/antiinstanton trajectories between the 
two ground states that are degenerate at ip = 7r and the parameter Y may be 
interpreted as the instanton fugacity. 

A dimensionless running coupling Y[ nst (L) is defined as 

Y lnst (L)=YL 1 -° . (47) 



From Eq.( T4T|) . one sees that, for g < 1, (i.e., when the Neumann fixed point is 
stable), instantons are a relevant perturbation, since Y[ nst (L) scales as ( -^-J 
However, the scaling ceases to be valid as Yi nst (L) ~ 1. 

For g > 1, using the boundary perturbation operator introduced in Eq. (j4"5l) . it 
is straightforward to compute the leading corrections to the Josephson current. 
For instance, when ip ~ 0, one may approximate the partition function in 
Eq.flU]) as 



cxp 



(Po,{0}\ 1 + 



2Y 2 J dr 1 dT 2 T r [:cos(^e(^ 1 ))::cos(-^e(ir 2 )):]}|P 0) {0}) . (48 



y/2 



By inserting in Eq. fT48|) the explicit form of the correlator 



(P ,{0}|TV[: cos 



cosh 



2 sinh 



2g 



(49) 



one gets the partition function Z[ 



cxp 



rWi-<f 



(3L 



Y 2 T[1 - 2g] 



r[g(i - £)] 



r[i 



r[g(i + £)] 
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Computing the logarithmic derivative of Eq. (|50|) with respect to <p and dividing 
it by —ft, one obtains 



e>pi y _ e , 2LT[l-2g] y2 g (■ r[g(l-f)] + T[g(l + £)] 



2ttL 
2ttL 



- ^ CM 



+ 



2nv pl d<p\r[l-g(l + £)] r[l - ^(1 - £)] 



(51) 



with 



G[g] 



2LT[l-2g] <9 2 r rfo(l-£)] 



+ 



4*%, a^ir[i-0(i + £)] r[i -<,(!-£)] 



.(52) 



For small 99, y, ifg provides a small correction to the slope of the Josephson 
current. As we shall see in detail in the following subsection, for ip = tt and 
g > 1, instantons at the Dirichlet fixed point remove the degeneracy between 
the minima and smooth the edges in the sawtooth-like Josephson current. For 
g < 1, of course, instantons are irrelevant operators. 



2.4 Instanton constribution to the Josephson current at the Dirichlet fixed 
point for ip = tt. 



In the previous section, it has been shown that, for ip = n, the state Pq, {0}) 
becomes degenerate with Pi, {0}). This corresponds to a discontinuous jump 
of e*gv p \/L in the Josephson current, between the values corresponding to 
the states that are degenerate at <p — tt. Since instantons connecting the two 
states should remove the degeneracy, for Y 7^ one should expect that the 
Josephson current becomes a continuous function of (p at (p — tt. 

In order to compute the F-dependent corrections to the ground state energy, it 
is useful to partition the Hilbert space into subspaces on which is diagonal. 
If \Pk, denotes a generic state in the sector corresponding to the eigenvalue 
Pk of the zero-mode operator, one defines the states belonging to the above 
subspaces as 

\±,*) = -^[\P ,*)±\P 1 ,V)} . 
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Thus, the partition function at tp = ir may be approximated as 

P P 
Z d [tt] =Tr|_ )$> |e-^T T exp[| dr # B (zr)]} + Tv ]+>9) ^ Hd T t exp[| dr H B (ir)]} 



where Tri T $\ denotes tracing over the sector of the Fock space spanned by 
the states |±, 

At low temperature, Eq. fl53l) may be approximated as 

P 



Z D [n] » (-,{0}|e-^T r exp[| dr H B (ir)]\-, {0}) 



o 

P 



+ (+,{0}|e-^T T exp[| dr H B (ir)]\+, {0}) , (54) 

o 



where |±, 0) are, of course, the ground states of both subspaces. By expanding 
T r exp[Jo dr H B (ir)] in a power series of Y, one gets 



P 

(±, {0}|e-^T T exp[ J dr H B (ir)}\±, {0}> 



o 

P 

\ ,.'/7n'..l" ^ 1 

exp 



R 9™y\ 

P 4L 



OO 1 M M 

E -T / rfr i--- / rfr m (±,{0}|T r [ J ff B (ir 1 )... J ff B (ir m )]|±,{0}) .(55) 
m=o 



When computing the partition function, one has to sum over contributions 
from connected, as well as from disconnected diagrams. By applying Wick's 
theorem, one obtains 

M 

(±, {0}\T T [H B (i Tl ) . . . H B (ir m )}\±, {0}> = n<±, {0}\H B (in)\±, {0}) + 

i=i 

M 

Y: (±,{0}\T T [H B (tr n )H B (tr n )]\±,{0}) ]J (±, {0}\H B (i Ti )\±, {0}> + . . . .(56) 

ii 7^2=1 i¥=3i,h 



Since (±, GS|.Hb(zt)|±, GS) = ±F, it is straightforward to evaluate, using 
Eqs.(155"l 1561). the leading F-dependent contribution to the energy of the states 



19 



±, {0}); namely E± is given by 

_ gnv pl L f T[f] r_[fg] | 

4L TT^pi " J if[T3|^ + f[r^fii ' (57) 



and, thus, even to the first order in Y, instantons remove the degeneracy 
between the ground states. 

The observable consequence of the removal of the ground state degeneracy is 
the smoothing down of the Josephson current, which becomes a continuous 
function of the applied flux <p when it crosses the value (p — tt. 

In order to see this effect, one has to resum over instanton contributions to 
the partition function, by setting ip = n + e, with e/7r <C 1. In order to obtain 
a meaningful result, one has first to sum over the instanton contributions, and 
then to compute the leading corrections in e. 

When ip = 7r + e, the low temperature partition function may be approximated 

as 

Z D fr + e] = ]T Tri^Je-^T.expf-c^e]]) , (58) 

a=± > 



with 



SS[Y,e] 



ft 

2Y J dr cos 
o 



'e(ir) 



gvpi 



eP 



9£_ 
4ttL 



(59) 



When e = 0, one obtains the result of the previous section. When computing 
the leading corrections in e, one has to notice that P mixes the two subspaces 
labeled by ±; namely 

(±,*|P|=f,*> = =F^ • (60) 



Thus, the leading e-dependent contribution to Eq. (jo9j) . d>Z£>[e], appears only 
to the order e 2 , and is given by 



^ r i ^ r i -fl^Ei f 2 \ - ^ f (m + l)(m + 2) / gv p ie\ 



2 

X 
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(a,{0}|(/3P) 2 |a,{0}>n / d n (a, {0}\H B (irj)\a, {0}) 

- s»w^ £{ i =^i^(w-[i + (-in 



XI eX P 



a=± 



2 A 



(61) 



To the leading order in e, the ground state energy derived from Eq. flBTj) is 
given by 

Eo ^ Eos ^ 0] -y-(^§ + ^ . (62) 



Eq. (1621) may be regarded as the lowest order term in the expansion of the 



function — yY 2 + (^r) • Although our approach allows to compute, in prin- 
ciple, also higher-order corrections, for our purposes it is enough to notice 
that, already to the leading order in e, one gets that the current behaves as 
a continuous straight line, with no sharp, discontinous, jumps, at e = (i.e., 

(fi = 7r). 

* dE GS [e] gv pl f gv pl \ 2 e 



Eq. fl63l) is not only consistent with the results of previous analyses of rf- 
SQUIDS [7f51] . but also with a strong coupling expansion of the ground state 
energy of the two-boundary sine Gordon model derived in Ref.JTJ] using the 
Bethe ansatz approach, once one of the two couplings is set to oo, and a strong 
coupling expansion in inverse powers of the other coupling has been performed. 
This makes us confident that the smoothing effect of the Josephson current 
due to instantons -exactly computed in Ref. [15] for g — 1- is effective for any 
value of the Luttinger parameter g. 



3 The dc-SQUID 



A dc-SQUID is realized as a superconducting loop, interrupted by two Joseph- 
son junctions, of nominal values Ej t i and Ej y2 - It is possible to generalize to 
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4> 



L 




4> 



R 



Fig. 1. Sketch of the dc-SQUID device with the leads. 

the dc-SQUID the renormalization group analysis developed so far for the rf- 
SQUID, obtaining that, also in this case, the only alleged fixed points are a 
weakly coupled (Neumann), and a strongly coupled (Dirichlet) one. 

The device is drawn in Fig. (CD): it is made out of two superconducting wires, 
connected by two Josephson junctions, of nominal values Ej^Ej^., with a 
magnetic flux $ piercing the loop. $2(2;) are the plasmon fields in the 

arms connecting the two junctions. 

When, in the following section, we investigate the transport properties of a 
dc-SQUID, the leads to which the device drawn in Fig. ([I]) is connected, are 
described by two additional plasmon fields and 

The Hamiltonian describing the dynamics of $1(2?), $2(2;) is given by 



H dcfi - — 

An j=l,2_ L 



dx 



1 



~~dt 



dx 



(64) 



while the Josephson energy of the two junctions is described now by the bound- 
ary interaction Hamiltonian H^ j given by 



H dc j = —Eji cos 



$ 2 (L)-$i(L) + | 



- Et 2 cos 



$ 2 (-L)-$ 1 (-L)-| 



with ip = $/$q. 

Again, it is most convenient to introduce the even- and odd-fields, $j >e / (x) 
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as 



^j,e/o\X) 



V2 



[$j{x) ± ; < x < L 



(66) 



By definition, the $j je (x)'s obey to Neumann boundary conditions at x = 0, 
while the $j )0 (x)'s obey to Dirichlet boundary conditions; namely, one has 



d$l, e (0,t) a$ 2 , e (0,t) 



<9x 



<9x 



; $i, o (0,t) = $ 2 , o (0,t) = . 



(67) 



The boundary Hamiltonian Hd c ,j gets contributions only from two linear com- 
binations of the fields <3> e , $ . Defining 

X(x) = -^[$ 2 , e (x)-$ lie (x)] ;((x) = -j=[^(x)-$ li0 (x)] , (68) 



it is easy to realize that the device is described by the reduced Hamiltonian 
H dc given by 



H dr = — / dx 



2n 



i foxy fdx s 



v p i \dt / 



— Eji cos 



- Sjo COS 



dx J 
(69) 



if^c in Eq. (l6"9"j) may be regarded as the two-field generalization of the two 
boundary sine Gordon Hamiltonian describing the properties of an rf-SQUID. 
By construction, at the inner boundary, x — 0, X(x) obeys to Neumann 
boundary conditions, 8X q^ = 0, while £(x) obeys to Dirichlet boundary con- 
ditions, £(0, £) = 0. The boundary conditions at the outer boundary, x = L, 
are obtained by requiring that the energy functional is conserved; namely, by 



gv pl dX(L,t) 



7T 



dx 



+ Eji sin 



X(L,t)+£(L,t) + | 



+ Ejo sin 



X(L,t)-£{L,t) 



gv p id£(L,t) 



dx 



Eja sin 



X(L,t)+£{L,t) + £ 



Ejo sin 



X(L,t)-£(L,t) 



For Ej y i = Ej y 2 = 0, X and £ obey to Neumann boundary conditions at x = L: 
= = 0. For -Ej,!, E J2 — >■ oo, the fields obey to Dirichlet boundary 
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conditions, t) +f (L, t) - f = (mod 2tt), X(L, t) -£{L, t) + § = (mod 
2tt). 

In the following, we shall use H^c to derive the renormalization of the Joseph- 
son energies and the functional form of the Josephson current in both the 
weakly coupled and the strongly coupled regimes accessible to a dc-SQUID. 



3. 1 The weakly coupled fixed point 



At the fixed point Ej^ = Ej$ = 0, X(x) obeys to Neumann boundary con- 
ditions at both boundaries, while £(x) satisfies mixed boundary conditions 
(namely, Dirichlet boundary conditions at x = and Neumann boundary 
conditions at x — L). The mode expansion of the field X(x,t) is, then, given 
by 

X(x, t) = q x + — + i\ - Y] cos(k n x) exp(-ik n v p \t) , (71) 

L 9 n 



with [ax(n), cex( n ')] — n ^n+n',o] Qx is the constant zero-mode operator and 
k n = jn. Px is, of course, the "dual" zero-mode operator, whose eigenvalues 
are the winding numbers. 

The mixed boundary conditions for the field £,(x,t) are satisfied by setting 

£{x,t) = -^Y^^sm[K n x}e tK ^ t , (72) 
V 2,7^0 n 



with [ag(n), c^(n')] = nd n+n i fi} and if„ = ~ (n + §J. 

By taking into account Eqs.( i7TT72l) . the boundary interaction may be normal 
ordered as in section 2. In particular, one obtains 



H dc ,j = -E JA : cos 



X(L)+£(L) + 



: -E J2 : cos 



X(L)-t(L) 



with E Jd =(^)'Ejj,j = 1,2. 

Again, the scaling equations for the pertinent running couplings are obtained 
from the requirement that the partition function Z is independent of the cutoff. 
Z is given by 

Z = / DID({e^ T T [e" Ss ]} , (74) 
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with 
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Sq = — / dr / dx 

2,71 J J 




1 fdX(x,iT)\ fdX(x,ir) 



dx 



9 



+ — / dr dx 
2tt J J 





v pl \ dr 



d£(x, it) 
dx 



(75) 



and 



S 



B 



dr [ e ^V hl (tr) + e-^V^iir)] 



E 



.1.2 



dr e^[y 1 ,_ 1 (zr) + e-^y_ lil (iT)] 



(76) 



while the vertex operators j,(t) are now given by 

Vafiiir) =: exp jz [&X(zr) + o£(it)] - a^j : ; a, 6 = ±1 . (77) 



To determine the scaling dimensions of the relevant operators, one needs to 
compute the correlators of the vertices written in Eq. (T77j) . They can be com- 
puted as discussed in appendix B; however, for the purpose of deriving the 
perturbative RG equations, one only needs the short- distance limit of the 
vertex correlators given by 



(TV [V a ,b (ir) V a > ,b> (ir')}) 



7™pli 

'r — r 



L 



(aa. — 1) 



<Wb',0 • 



(78) 



Eq. flTSl shows that the scaling dimension of V a ,b is h a ,b = l/g- 

Since £(r) does not have a zero mode, the O.P.E.'s generate, just as it happens 
for an rf-SQUID, higher periodicity terms, according to 

-h a , b -h a , b ,+h a+a , h+b , 

V a+a , >b+b ,{iT') .(79) 



V a ,b(ir)V a , 



[IT 
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To the second order in the couplings Ej t i, Ej t2 , Eq. (!79l) requires to add to Sb 
an extra counterterm given by 



5S B = -A 2 J dr [V 2>0 (ir) + V- 2 ,o(ir)] . 



JO) 



The running coupling constants Gi,G 2 are now defined as 



I _ 1 

L \ a aEjo 



(81) 



Since the vertices V± 2 ${iT) = e ±lip : exp [±2i£(r)] : are less relevant than the 
V±i,±i's, they do not contribute to the evaluation of the leading terms of the 
perturbative j3 functions for Gi,G 2 , 



dlnGi(L) 
d\n(L/L ) 



Pi(Gi, G 2 ) 



Gi(L) 



and 



d\nG 2 (L) 
d\n(L/L ) 



G 2 {L) 



(82) 



As for the rf-SQUID, one finds that both couplings are irrelevant for g < 1; 
namely, the Neumann fixed point is infrared stable, and the theory is pertur- 
bative in the couplings Ej t i, Ej t2 . 

When g > 1, instead, both couplings Gj(L) increase, when L increases, as 
(L/L ) 1_ ». Thus, just as for the rf-SQUID, the theory becomes nonperturba- 

tive as soon as L ~ L*, with L* ~ 2na [—^ — ) 9-1 , being Ej m i n the smaller 
coupling, between Ej t i and Ej >2 - 

To evaluate the Josephson current for g < 1, we compute the partition function 
for finite (p. The result is 



Z[<p] = Z (T T exp 



dr 



Ej i cos 



X(ir)+C{ir) + 



V 3 



+ E J2 cos 



X{ir)-C{ir) 



») 
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with 



n=0 



2g j irv ri ] „ 

-/3 pl m 2 



and q = e 13 l . 

Since the field -X"(ir) has a constant zero- mode term, a perturbative expan- 
sion in the boundary interaction yields nonzero contributions to the partition 
function starting only from the second order in Ej^,Ej t2 , implying that, at 
low temperatures, the leading 92-dependence is given by 



Z[if] = Z ({l + ~ j dn J dr 2 T T [H dC; j(n)H dC! j(T2)}]) 



p 

Z ({l + /3£j,i£j, 2 cos(^) J dr 



2 cosh ( ff t) 



2L - 

ITU 



J5) 



with 



B 



2§r[i + 



i6) 



result, the Josephson current is given by 



J $ = h m ^ 



2L 

e* — BE JA E Jj2 sm((p) 

7TU 



J7) 



The functional form of /[$] is similar to Eq. (|24p . but it is now proportional to 
Ej^Ej^. For g < 1, Eq. flSTl) is valid for any value of L, since the quantum cor- 
rections are perturbative in Ej^Ej^- When g > 1, instead, higher harmonics 
operators become more and more relevant as L increases, inducing a crossover 
to a sawtooth-like functional form of the Josephson current, as L approaches 
L* ■ 
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3.2 The strongly coupled fixed point 



In this subsection the strong coupling (Ejj/ (2ttv p i/ L) <C 1) regime of a dc- 
SQUID is analyzed and the explicit functional forms of the partition function 
and of the Josephson current are derived. To compute the partition func- 
tion with Dirichlet-like boundary conditions at the outer boundary, X(L,t) + 
£(L,t)+§ = (mod27r)), and X(L, t) -f(L, t) - f = (mod2vr)), one notices 
that 



X{L,t)=n{n 1 + n 2 ) ; £(L,t) 



n(ni-n 2 ) ■ 



with ni,Ti2 relative integers. 

Eqs. pSjl . together with the boundary conditions at the inner boundary, lead 
to 



2tt 



P^x 



— ^2 s i n 

9 n^O 



-x 



n 



(89) 



with [a^(n), a^(n')] = n5 n+n ' t o, k n = nn/L] the 0-mode operator of the field £, 
P%, has the eigenvalues 



(p m 
4^ + IF. 



(m = n\ — n 2 ) 



(90) 



Similarly, one obtains 



X(x,t) =q x + i 




n + 



x 



ut 



(91) 



with [a x (n),a x (n')\ = n8 n+n > }0 , K n = | (n + ±J , q x = J\ 7r(ni + n 2 ). 

As for the rf-SQUID, the strongly coupled fixed point is described in terms 
of the (now two!) dual fields Q x {x, t) and Q^(x, t), whose mode expansion is 
given by 



e x (x,t) = -2nP x -2^g-J2 



sin 



neZ 



7T / 1 

l{ n+ 2 ]X 



ax(^) e -ij(n+i)«* 



n + 



i 



,(92) 
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where Px is the canonical conjugate operator of qx, and 
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Bt(x, t) = 9 + —gP^ut + i2y/g ^ cos 



Tin 



-x 



a i\ n , 
n 



e l 



(93) 



Following the same steps used in subsection (2.2), one has that 

H D [Q X , 9 € ] = ^^(P f ) 2 + Y,M-n)a x (n - 1) + a^-n)a^n)} , (94) 



from which it is straightforward to compute the partition function at the 
Dirichlet fixed point: 



n=0 



1 



1 



1 _ „n+\ Ml- q n+l 



$>xp 

J kez 



Eq. (l95l) shows that, for if ^ 2kir + ir, with k integer, the ground state is 
characterized by a nonzero value of Pg, proportional to [ip\. Tunneling events 
from this ground state to the nearest (in energy) ground states are suppressed 
by the exponential factor exp[— (3irv p \g/2L]. 



Taking the logarithmic derivative of Eq. (l95l) one easily derives the Josephson 
current, which is given by 



rr i * r 1 d\nZ D [(p] e*gv p i 
I[cp] = -e Jim — = -^f^W ~ m) 



/3^oo (3 dip 



(96) 



From Eq. (|95p . one sees that the degeneracy among the minima of the Joseph- 
son energy is removed by the magnetic energy E$ = 29 ™ pl (— ^ + • Since 

E$ = E^ +1 , for <p = 7t + 2vrm also a dc-SQUID may operate qubit be- 
tween these two levels. When the couplings are large, but finite, phase slips 
in the plasmon field (instantons) will induce a renormalization of the physical 
parameters, as we shall show in the following subsection. 



3.3 Instantons and boundary perturbations at the Dirichlet fixed point 



In order to construct the leading boundary perturbations at the strongly cou- 
pled fixed point, one may use again the DEBC-method [19] . applied to the 



29 



primary bulk operators involving the dual fields Qx, ©£■ A generic, imaginary 
time, bulk operator at the strongly coupled fixed point may be written as 
Va,p(i T ) =: exp{z[a0x(«T) + pO^ir)}}, where the imaginary time boundary 
fields are given by 

®x{it) = e x (L, it) = -2nP x - 2^g £ e -£(«+*)- , (97) 

ndZ 71 + 2 



and by 

e^T) = e € 0MT) = e + ^2gP^uT + t2^Y,(- 1 T—e~ TuT -( 98 ) 

J— j / n Tl 



From the commutation rules between the zero modes of the two fields, one 
obtains 

[qx,V a , p {iT)\ = 2™V a ,p{iT) ; [P 6 V Q)/3 (zr)] = 2ir/3V atP {ir) . (99) 



Thus, V a fi(ir) changes the eigenvalue of qx by 2na and the eigenvalue of 
by 2nj3. Using Eqs. ()97|l98|l9"9~j) . one may write the leading boundary operators 
as 

V a , b {iT)=:exp{~[ae x {iT) + bez{iT)}\ ; a, b = ±1 , (100) 
Vx )± (zt) =: exp{±iQ x {ir)} : ; V^,±(ir) =: exp{±i6 € (ir)} . (101) 



Vi,i(ir) and VLi -i(ir) (Vi,-i(*Y) and VLi^ir)) change by ±27r the phase at 
junction 1 (2), while they leave the phase at junction 2 (1) unchanged and 
thus they describe phase slip operators at junction 1 (2). Similarly, Vx,±(ir) 
(V^±(iT)) change the eigenvalue of the overall (relative) phase X(L) (£(£)) by 
±2n leaving the eigenvalue of £(L) (X(L)) unchanged: thus, Vx,±(ir) do not 
affect the Josephson current across the ring, while V^±(ir) provide corrections 
to the Josephson current when if = 2kn + n. In the following we shall refer to 
Vx,±{ir) (V^±(ir)) as q- (it- ) Instanton operators. 

The leading boundary perturbation at the strongly coupled fixed point is given 
by an arbitrary linear combinations of the operators listed above, as 

H dCl j = YtiVt^ir) + V-i,-i(iT)} + Y 2 [V 1 ^ 1 (ir) + VLi,i(zr)] 
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+ Y x \V x>+ (ir) : +V x ^{ir)} +Y^ + {ir) : +V { ,_(zr)] 



(102) 



From the mode expansion of the dual fields, it is straightforwad to compute the 
(low-temperature) correlators among the boundary vertex operators appearing 
in Eq. fll02p . For \(p\ < n, for instance, one has that 



(T T [7i,i(ir)y_i_ 1 (ir')]) = (T^y-i^ir)^,-^/)]) 



e 4l < > 



sinh 



29 ' 



(103) 



(T T [V +;X (ir)V-, x (ir')}) 



cosh 



sinh 



19 



(104) 



and 



{T T [V + ,Sr)V^{ir')\) 



2 sinh 



2L 



4fl 



T — T 



(105) 



Eqs.f ll03lll04|105p imply that the fugacity of the q- and 7r-instantons and of 
the phase slips at the junctions scale with the SQUID size as 



Yx(z){L) = Y m) (v p iT x ) 



L 



l-2g 



v p \T x 



Yi(2)(L) = Yi(2)(v p iT x 




.(106) 



The "instanton size" T x in Eq. (11 06 ft is estimated to be of the order of J gL/[irv p i(Eji + Ej 2 )] 

m 

Eq. (ll06p shows that, for g < 1, phase slips at the two junctions provide the 
most relevant perturbation to the Dirichlet fixed point, rendering the strongly 
coupled fixed point not infrared stable. For finite (3, L and for <p 7^ 2rck + 
ii, however, the phase slips are always suppressed by the magnetic energy 
associated to the eigenvalue of P^. At variance, when g > 1, 7r-instantons 
are an irrelevant perturbation which, for L < L^L^i (i.e., of the healing 
lengths associated to the two junctions), only smooths down the edges of 
the sawtooth shape of the Josephson current: for L — > 00, the ground state 
exhibits a discrete Z-symmetry. 
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4 Transport in SQUID devices 



To probe the different regimes attainable by SQUID devices, it is most useful 
to look at their dc-conductance. To get conduction of current across a SQUID 
device, one should be able to connect it to two leads, enabling to apply a 
biasing voltage V. This can be achieved easily for a dc-SQUID while, for the 
purposes of this section, it is most convenient to regard the rf-SQUID as an 
inhomogeneous chain of Josephson junctions connected to two bulk super- 
conducting leads at fixed phase difference cp and at finite biasing voltage V; 
the chain is made by junctions of equal strength with a weak link of nomi- 
nal strength Ej located at its center. For V = 0, this inhomogeneous chain 
mimicks the response of an rf-SQUID to an external magnetic flux (p [21]. In 
this section we study the trasport properties of both devices, evidencing the 
different current (normal and Josephson) response at weak and at strong cou- 
pling. Our boundary field theory approach well reproduces the known results 
of Refs. [Tf5] and, for a dc-SQUID, evidences interesting interference effects 
between the current flowing through the two junctions. 



4-1 The Josephson junction chain at the weakly coupled fixed point 



The relevant dynamics of a Josephson junction chain with a weak link con- 
nected to two bulk superconductors at fixed phase difference (p has been al- 
tready analyzed in Refs. [3~5f2T] . The properties of this device may be described 
by the bosonic field = — x)], where $>r(x) and are 

the fields describing the phase of each junction in the right and the left half- 
chain, respectively. The current operator across the junction is then given by 
j(t) = — 9e ^ 1 d ^Q^ ■ Energy conservation provides the dynamical boundary 
conditions for Q(x, t) at the outer boundary, which is given by 

V ^^^ + V2Ejsm[V2^L,t)+<p] = ; (107) 



using Eq.f ll07p . the current operator may be written as 

j(t) =e*E J sin[v / 2$(t)+^] ; ($(*) = $(£,*)) . (108) 



At the weakly coupled fixed point (Ej = 0), $(x, t) satisfies the Neumann 
boundary conditions at the outer boundary; namely, = 0. 

In order to apply a voltage V to the device, one may bias at voltage V/2 and 
— V/2 each one of the two leads. In the small ^j-limit, it is safe to assume that 
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each chain is at equilibrium with its own lead, and that the two half chains 
are disconnected. Accordingly, the voltage bias V is introduced by adding to 
the total Hamiltonian the term 

e*Vg f , d$ 



H v = / dx — , (109) 

irv p i J at 



which can be accounted for by a mere shifting of the field $, as 

<$>{x, t) -»• ${x, t) - e*Vt . (110) 



As a result of Eq. (II 101) . the Josephson current acquires an explicit dependence 
on time and, to the first order in Ej, reads 

Ij[(p,t]=e*EjSm[ip + e*Vi\ . (Ill) 



Thus, in addition to the Josephson current, there is a normal time-independent 
dc-current which, as it will be shown in the following, becomes the leading con- 
tribution for a pertinent choice of the infrared cutoff leading to the celebrated 
Kane-Fisher formula [9]. 

Since the normal current is independent of time, it may be computed using 
the imaginary time formalism. From Eq. flllOl) . one obtains 

j (it) = e*Ej : sin[V2$(zr) + <p + e*Vir] : , (112) 

with 

$(ir) = -^Et" 1 )"^- e-T( n HW- . (H3) 



n , 2 



From Eqs. flll2|113p . it follows that the leading contribution to the normal 
current appears only at the second order in Ej, since 



I w j.{Ejf J dr' (T T [: exp[V2$(ir) + if + ie*Vr] : : exp[-V2$(ir') - <p - ie*Vr] :]) 
o 

-j-^Ej) 2 J dr' (T T [: exp[-V2$(ir) - cp - ie*Vr] : : exp[V2$(i<r') + <p + ie*Vr'\ :]) 
o 
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2L 



TCV- 



pl 



dx 



sinh [ax] 
2 9 sinh f (x) 



2e*7i(Ej) 2 L 

TTVpi 



(114) 



with 



2Le*V 

TTVpi 



(115) 



and 1 36 1 



Ifal 



dx 



sinh [ax] 
[2 sinh(x)] b 



r[i 



r[-f + i] 



r[i 



r[§ + |] 
r[i + §-i] 



.(116) 



From Eq. flll4p . one sees that the dc-current depends crucially on the ratio a 
between the energy window defined by the applied voltage, e*V, and the level 
spacing characteristic of the device, nv p i/(2L). In fact, if e*V < ttv p \/L, one 
may expand Eq. (I114[) to first order in a, getting 



e*(Ej) 2 LT[l-^ 



2v 



pi 



(117) 



with 



B 



d_(Tl 
dz\T[l 



- + -} 

2 Z 9 l 



£[f+j] 
r[i + f-J] 
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2 = 



which yields the (normal) dc-linear response, to an applied voltage V, of a 

junction with renormalized strength Ej = (^jf^j 3 ] this is expected, since, in 
this limit, ttv p \/L provides the infrared cutoff. If, instead, v p i/L <^ e*V, one 
may compute the current in a large-a expansion, using the Stirling formula, 
T[z] ~ ^/2ne~ z z z ~2 , which yields 



e\Ejf L 

7TV p i 



sin 



71 


r 






'2Le*V~ 








g_ 




9_ 




TTVpl 



±-1 

3 



(119) 



Eq. flll9p is the celebrated Kane-Fisher formula for the current across a con- 
striction in a spinless Luttinger liquid. Our result is, after all, not surpris- 
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ing, since the wires connecting the junctions have been regarded as one- 
dimensional spinless Luttinger liquids, and the junctions as boundary interac- 
tions. It should be noticed, however, that the normal current is proportional 
to (e*V)s since e*V is the pertinent infrared cutoff when an external bias 
voltage is applied to the device. At variance, in the computation of the Joseph- 
son current, the pertinent cutoff is provided by the plasmon's level spacing 
s = 2-KVpi/L. Thus, in the limit where s <C e*V, our analysis shows that the 
leading contribution to the total current across the junction is given by the 
normal current. 



4-2 The Josephson junction chain at the strongly coupled fixed point 



In this section, the current across the Josephson junction chain is computed 
at zero phase difference tp as Ej/(ttu/L) — > oo. Here, it is most convenient 
to resort to the dual formulation, based on the dual field Q(x,t). The total 
Hamiltonian (including the leading boundary perturbation) is given by 



H\ 



An{2g) 



L 

Jdx 




+ u 




Y[: e 7f e(L) : + : e ^ e(L) : ] , (120) 



where the dual field Q(x, t) obeys to Neumann boundary conditions at x = L. 
Thus, the mode expansion of the field 0(ir) = <d(L,ir) is given by 

e(ir) = 0(L,ir) = 6 + i^2gPur + i2^^(-l) n ^ e - s ? UT . (121) 



Since at the strongly coupled fixed point, the two half chains cannot be re- 
garded as isolated from each other, one must apply a voltage V by biasing 
the right-handed chiral mode of <$>(x,t), with respect to the left-handed ones. 
This corresponds to adding to H[Q] a volt age- dependent term given by 

e*V ) , 9$ e*V f , dO „ nn . 

p 



The imaginary time current operator is 
5S E [Q;V,a] 



Sa(r) 



a=0 



(123) 
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where Se[®', V, a] is the Euclidean action with a source term for the current, 
given by 



S E [&;V,a] 



An(2g) 



L 

dr / dx 



o o 



v p i \ dr 



dx 



2Y / dr : cos 



egr) 



. e*V r , / , dQix.ir) ie* 
— i I dr I dx 



2ttv 



pi 







<9r 



2ttL 



dr / (Wr) *"^ , (124) 
or 



o o 



By shifting the field according to 
Q(x,ir) Q(x,ir) 



s/2 \[2 



— 2ige Vt — % — —^-A[j) , 



(125) 



with A'(t) = a(r), one finds that the current is given by 



5S E [G;V,a] 



5a(r) 



a=0 7T 



-y + 4y- 



■ sin 



V2 



ie*Vr) : .(126) 



From Eq. (11261) one sees that, apart from the constant term g ^ V, the current 
may be computed as in the previous section, provided that one substitutes A/g 
with g and Ej with Y . Thus, to the second order in Y one obtains 



/ = £ (£)! v + 2 f z!£ sinMr[1 _ 29] 



7T 



7T 



2Le*V 

7lV p l 



2g-l 



(127) 



As it happens also at the weakly coupled fixed point, the celebrated power- law 
dependence of / on the applied voltage V emerges when 2Le*V/(iTVp\) 3> 1. 



4-3 The dc-SQUID at the weakly coupled fixed point 



For a dc-SQUID, it is most convenient to represent the leads as two quantum 
wires with the same Luttinger parameters g and v p \, with plasmon fields re- 
spectively given by $l and The two leads are assumed to be connected 
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to the dc-SQUID with smooth junctions, that is, three-wire junctions with no 
backscattering in any arm. Thus, one should describe the SQUID's arm with 
four fields, rather than with two; namely $ Uj l, & u ,r for the left-upper and for 
the right-upper arm of the ring, $d,L>$d,R f° r left-lower and for the right- 
lower arm of the ring, respectively. Since the connections between leads and 
SQUID's arm are assumed to be ideal, they are described by the strongly cou- 
pled fixed point of the three quantum wire "Y"-junction studied in Ref.|19j. 
with no concatenated flux. Thus, by centering the Y-junction at x — 0, the 
boundary conditions are given by 



d_ 

dx 



%(X, t) + $ u , a (z, t) + $d,a(x, £)] 







IE=0 



(128) 



which embodies the current conservation at the connection, and by 

- 2$ a (0, t) + $ Uji (0, t) + $ d>a (0, t) = $ U)O (0, t) - $ d , o (0, t) = , (129) 

with a = L,R. 

The relevant fields for the Josephson junction dynamics are given by the "odd- 
parity" combinations 

&u/d(x, t) = -j=[® u/d>R (x, t) - $ u/d , L (x, t)} ; ${x, t) = -^=[<& R (x, t) - $ L (x, t)\ , (130) 



which are related to the fields £(x,£) , X(x,t) introduced in section (3) by 

<f> u {x,t) = ±[X(x,t)+t{x,t)] ; $ d (x,t) = ^[X(-x,t)-Z(-x,t)] .(131) 



As one sees from Eqs. fU28|129j) . introducing the contacts does not affect at all 
the boundary condition of £(x, t) at the inner boundary (x = 0). The boundary 
conditions for X(x, t), instead, explicitly depend on the applied voltage bias, 
whose effect may be accounted for by including in the total Hamiltonian a 
pertinent "source" term. 



At the weakly coupled fixed point, biasing the left-hand lead at a dc- voltage 
V/2 and the right-hand lead at a voltage — V/2, corresponds to adding to the 
SQUID's Hamiltonian a voltage dependent term given by 



9 



V2e*V 



v 



4tt 



dx 



d$ R {x) d$ L {x) 



dt 



Ot 



ge*v 

2tt 



dx 



dX{x) 
dt 



(132) 
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which is equivalent to a shift of <$>(x,t), which is linear in time, namely, the 
field $(x, t) should be replaced by t) - e*Vt. From Eqs. (jl28|129|132p one 
gets that also $> u /d(x,t) must be shifted as 



e*V 

&u/d(x,t) -> $ u/d (x,t) - -—t 



(133) 



As a consequence, one sees that X(x,t) takes the additional contribution, 
given by 



X(x,t) -*X(x,t) -e*Vt 



(134) 



which introduces in the boundary interaction Hamiltonian an explicit depen- 
dence on r, given by 



H dc ,j -> H dCyJ (iT) 



-Ej i cos 



X{iT)+OiT) + ^+ie*Vr 



Ej o cos 



X{iT)-i{ir) - ^+ie*Vr 



(135) 



at finite V, the current operators across the two junctions are then given by 



j'i(ir) = e*E J1 ean 



X{ir)+^{ir) + ^ + ie*Vr 



(136) 



and 



j 2 (ir) = e*E J2 sin 



X(ir) - Z(ir) +-- ie*Vr 



(137) 



As in the previous section, the current across the two junctions is computed 
from a perturbative expansion in Ej t i, Ej^- One obtains the following results 

• Current across the junction 1: I\ 



h = h,i + I 



1.2 
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Ii,! is given by 

P 

h,i = e*(E JA ) 2 J dr' (T T {:sm{X(iT)+Z(iT) + ^ + ie*VT] : x 



: cos[X(*r') + £(ir') + | + ie*W] :}) = y (£j,i) 2 / dr { 

o 

which, for ALe*V/ \v p \ti) > 1 yields 



sinh [e*Vr] 
2sinh(^r 



2e*(E 7 , 1 ) 2 X . 
J 1 1 « ! sin 



7T 












_<7_ 







(138) 



J ll2 = e*Ej, x E J>2 I dr 1 (T r {: sin[X(2r) + £(ir) + | + eTr] : x 



: cos[X(zr') - £(zV) - | + e*Vr'\ :}} 



e*Er A E, 2 / dr 



3 e* Vr+iip 



cosh (^r) 
e i¥, J[a] - e -^J[- a 



cosh (^r) 



with 



J [a] 



COsh a [xj 

where a = 4e*VL/ (v p \7i) and 
J[o] = 2l~ 1 e- i <^)B 



(139) 



(140) 



/a 1\ 2 



(141) 



where the incomplete - Beta function is defined as B[z, a, (3] = Jq dtt a *(1 — 

For large a and 1 — 2/g > 0, one may consider the asymptotic expansion 
of the incomplete- Beta function [36] , 



-i;- 






;i- 


2" 




7T 


.7T 
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6- 






w z exp 


—i—a 


— z— 


(H 














2 


2 







5 

from which 

2e*Ej,E Ji2 L 



--u 



irv p i 



r[i 



COs(<£>) COS 



2e*1/L 7T 

+ - 
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2e*VL 7i 




~4Le*V~ 


+ - 


} 




v P i g 




7lV p l 



+ sin(<yj) sin 

Current across the junction 2: I 2 
h — h,i + h,2 1 



(143) 



The two contributions are computed below: 



I2.2 = < U-:.,.2) 2 I dr' (T T {: sm[X(tr) + £(zr) + | - ie*Vr) : x 



ry/. A , t/ . , ■ * T/ /] ~l \ 2&* {E J 2 ) 2 L . 

: cosl A ur ) + 4(zr ) + — — xe Vr J :}) ~ sin 

2 1 TTVpi 



r[i--] 



4Le*U 



_ 2 



.(144) 



(2) 



J 2)1 = e*E JA E J>2 J dr' (T T {: sm[X(ir) + £(zr) - | + e*Vr] : x 



^osfX^rO-e^'rO + l + eW] :}) 



1 1 < cosfo?) COS 

™pl 9 



2e*VL vt 
+ - 



+ sin(yj) sin 



2e*1/L vt 
+ - 


} 


"4Le*\/" 


^ P i 9 







(145) 



The net transport current across the dc-SQUID is obtained as the sum of the 
currents across the two junctions, that is 



I = h + h 



2e*L 



TIU 



+ Ej A Ej t2 cos(v?) cos 



L --]{[(% 


r- 


H fe) 2 


] sin 


7T 


9 I 






_9_ 


~2e*VL tt" 
+ - 
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■4Le*y 


2-1 

9 




v P i 9 











(146) 



Eg. (1146[) contains an "interference" contribution (proportional to cos(^)), rep- 
resenting the cumulative effect of a bias voltage applied to junction 1 on the 
current flowing across junction 2, and of a bias voltage applied to junction 2 
on the current flowing across junction 1. This is a remarkable result of the ap- 
plication of the boundary sine Gordon techniques to SQUID devices. Treating 
the plasmon modes as modes of the dynamical fields, allows for an explicit 
computation of the interference contributions to the dc-current. 
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The circulating current across the dc-SQUID, J, may be computed as the 
difference between the currents flowing across the two junctions, Ij = I\ — 1%. 
The result is 



4e*E J1 E J2 L . 2. . , 

J = ■ ■ — r 1 sin(<£>) sin 

Ti^'pi g 



2e*VL 7i 




~4Le*V~ 


+ - 


} 




v P \ g 




TVVpi 



(147) 



which, again, comes from pertinently taking into account correlation functions 
of the plasmon modes computed using the boundary sine Gordon theory. 



4-4 The dc-SQUID at the strongly coupled fixed point 



At the strongly coupled fixed point, as long as one neglects phase slips at 
the junctions, one should not expect any resistance across the junctions. As a 
consequence, the external voltage V is added by biasing the left chiral mode 
of the plasmon field at each lead with respect to the right chiral mode, with 
opposite biases in the two contacts. As discussed in section (7.2) for the chain, 
this is accounted for by shifting the field <&(x,t) according to 



e*V 

${x,t) — > ${x,t) X 
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Due to Eqs. fll28lll29[ I148j) . one has that, at finite V, also $ u /d(x,t) must be 
shifted as 



e*V 

$u/dOM) -> ®u/d(x,t) - —X 



Pi 



(149) 



Eqs. (11311) imply that 



X(x,t) -> X(x,t) 



e*V 

Vpl 



(150) 



which is equivalent to adding to the total Hamiltonian a voltage dependent 
contribution given by 



e*Vg f , dX(x) e*Vg ) , de x (x) 
H v = / dx „ ' = -- — - / dx 



dx 



2-KVpi 



dr 



(151) 
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The current operators across the two junctions are now written as 

= jx(iT)±jt(ir) . (152) 



h,2{ir) = -i- 



dQxiir) ± dQs(ir) 



dr 



dr 



An useful representation may be derived by adding to the total Euclidean 
action Se[&x, a x, o^] the source term given by 



5S% OUTCC 



2ixL 



I dT I 





dx 



. ,dQ x (ir) . .dG*{ir) 



dr 



dr 



, (153) 



and by taking the functional derivatives of Se [©x, 8$,flx, og] with respect to 
ax, d£, after shifting the dual fields according to 

Gx(x,iT)^e x (x,iT)-2i^p±Ax(T)-ige*VT , (154) 



and 



e € (x,<r)-e € (x,<r)-2<^-^^(r)+<^T , 



(155) 



with A^(r) = a x (r), ^(r) = a 5 (r). 

The Euclidean action in the dual representation is then given by 



S E = Y — !— [ dr [ 



L 




J dx 










' d<d a (x,ir) 
dx 



—2Yi J dr cos 
o 





ge*V - v p i 



4L 



r + i^^[Mr)+Ax(T)} 



2Y 2 J dr cos -(-&x(ir) + O^ir)) - i 



geV + v pl — 



r + i^P±[Mr) - A x (r) 



P L 
+ A^2q)l ^ I dX { 







2ge*V + 2e> pl 



a x (T 



+ 



o * a d T ) . v 



.(156) 
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As a result of Eg. (11561) . one finds that the current operators jx, j$ may be 
written as 



]x{n~) = V + At j- — sin 

7T L 



(9x(ir) + e^(ir))+i 



ge*V + v p i 



e*gv pl Y 2 . 
2t sin 



T 



AL. 
(157) 



r 



and by 



/ ■ \ e*0u„i e*gv pl Y 1 
j£{tT) = A _ T ip + At ^ — sm 



AnL 



A .e*guY 2 . 
+ At — - — sm 



-(0 x (zr) + 5 (zr))+z 



; (-0 x (zr) + 6 5 (zr))-z 



AL 



T 



, (158) 



T 



From Eqs. (ll57lll58l) . one sees that, to the lowest order in the phase slip con- 
tribution (Yi = Y 2 = 0), the stationary component of the circulating current 
reduces to the value of the Josephson current at zero voltage, ^^f, while the 

conduction current is given by 9 ^ e J V, just as for a single chain at the strongly 
coupled fixed point. 



Higher-order corrections, including the effects of phase slips, may be computed 
following the same steps used in the previous section. The result is 



X 



Ae*g{Y y 



71 



r[i-2d 



r[i-if 



+ 



Ae*g(Y 2 ) 2 



71 



ni-2g] 



9} n^ + a --g] 



T[l + a -±-g] T[l- a -±-g] 



(159) 



with 



2ge*VL wp 

7TU n 



(160) 



and by 
8jz 



^e*g(Y 1 
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ni-2g] 



n- a -r+g] 



r[i 



9] r[i + ¥-g] 
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4e*<7(Y 2 



7T 



r[l - 2g] 



m + g] T[- a -±+g} 



r[l + 2±-g] T[l-^-g] 



(161) 



For if ~ and for \a±\ <C 1, one may expand Eqs. fll59|ll61l) to the first order in 
a±. In the symmetric case {Y\ = Y 2 ), this provides a correction to the slope of 
both the dc-current and of the stationary component of the Josephson current, 
so that one gets 

Sjxx + 2,1 W , (162) 

Z7T I 71U J 



and 



s ^m^ ij ^^-^ ■ <««> 



with 



^lr[i-z + ^] r[i + ^-^]' 



Eqs.( ll62lll63p show in a rather simple context the way in which relevant per- 
turbation may affect both the Josephson and the normal currents across the 
device. 



5 Comments and conclusions 



We provided a framework where two-boundary sine Gordon models may be 
used to investigate the relevant nonperturbative features, occurring in the 
renormalization of the Josephson couplings Ej ti of a multi-junction SQUID. As 
pointed out in Ref . [7] , this is needed since -for large fluctuations of the phase 
of the order parameter- the effective boundary conditions on the plasmon 
fields depend on the Josephson energy of the junction. Our analysis shows 
that, when g > 1 and L > L*, the Josephson couplings affect the dynamics 
of the plasmon field only as a boundary term, while, for L < they can be 
regarded as quantum impurities in a superconducting loop. For g < 1, instead, 
the Josephson couplings have to be always regarded as quantum impurities. 
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The boundary field theory approach to SQUID devices proposed in this paper 
turns out to be very powerful in providing a systematic procedure for com- 
puting not only the Josephson currents (and the leading corrections induced 
by the relevant perturbations), but also the dc-currents flowing in a SQUID 
device, due to externally applied voltages. For pertinent choices of the infrared 
cutoff, the results of the boundary field theory well reproduce what has been 
obtained in Ref . [9|7] . allowing, in addition, to appreciate new and closer con- 
nections between the theories of Josephson superconducting devices, quantum 
brownian motion, fermionic quantum wires, and quantum impurity problems. 
Furthermore, for dc-SQUID devices, our analysis points out the existence of 
new remarkable interference effects between the currents flowing through each 
junction. 



We investigated here SQUID-devices with one or two junctions. For these 
systems, we have found two stable regimes characterized by stable fixed points 
of the pertinent RG equations, driven by the strength of the bare Josephson 
energy. It would be interesting to ascertain if superconducting loops with 
more than two junctions may exhibit new (in addition to the by now well 
known weakly and strongly coupled fixed points) renormalization group fixed 
points at a finite Josephson coupling. Furthermore, it could be instructive to 
investigate the properties of a multi-junction SQUID device, for which g — 1, 
since the boundary sine-Gordon model provides an exact solution for this value 
of the Luttinger parameter. We feel that mult i-j unction superconducting loops 
may become an interesting laboratory for testing the physical properties of 
field theories describing interacting Luttinger liquids with pertinent boundary 
interactions. 
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A Scaling properties of boundary interaction strengths 



In this appendix we shall derive the scaling laws for boundary interaction 
operators. 

Let us consider a one- dimensional system described by the Euclidean action 

M P 

S = S -J2*i I dr^(r) , (A.l) 
1=1 i 



where So, the action for K independent massless Klein-Gordon fields, {4>j}, 
with pertinent boundary conditions, is given by: 



K 



So 



E / dr 



dx 



1 



dr 



+ v 



pi 



dr 



(A.2) 



The boundary operators {$^(r)} are functionals of the fields {<ftj} with scaling 
dimension hg. As a consequence, for r — > r', one gets 

Gt (r - O = (Mt)Mt'))° V _ X , (A.3) 



where (. . .)o denotes averaging with respect to the "free" action So, while L 
is the "size" of the system. 

It is well known that, in order to introduce scale invariant interaction terms, 
one needs to define dimensionless coupling constants, given by gt(a) = A^a 1- ^, 
where Xe is the coupling constant appearing in Eq. flA.ll) . and a is a short 
distance cutoff. To first order in the coupling strengths, the renormalization 
group equations are given by 

= ■ (A.4) 

dm{a/a ) 



Higher order corrections to the renormalization group equations come from 
nontrivial short- distance fusion rules of the 5>/s, which are represented by the 
O.P.E.'s 

MW) E W _%U-^ T) ■ (A - 5) 
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If the O.P.E.'s coefficients Ctv p are different than zero, the scaling equations 
get contributions which are of second order in the couplings. To derive them, 
one should start from the partition function 



exp 



M 



dr a' 



-l+h. 



(A.6) 



Upon introducing the cutoffs a and L and by expanding Eq. (IA.6I) up to the 
third-order in the running coupling strengths, one gets 



M L/v pl L/v pl 



^^1+ gm , f dr f dT'a- 1+h 'a- 1+h "9(T-T'-—)(^(r)^(r')) 



- L/v p i L/v pl 



M u "I "pi "pi 

+ 9i9i'9e" J dr J dr ' J dl ~" a 



II „-l+h( a -l+h e , a -l+h e i, x 







0{r ~r'- —)6{r> - t" - — )($,( r )$,,( r >,„( r "))o • (A.7) 

Vpl V p i 



If one rescales the cutoff, a — > (1 + e)a, with e< 1, additional contributions 
to the coupling constant renormalization arise, as a result of configurations 
where the arguments of two fields lie between a and (1 + e)a. The ensuing 
renormalization is derived by means of the identity 

9(t - t' - — (1 + e)) 9(t' - r" - — (1 + e)) « 9(r - r' - — ) #(r' - r" - — 

V v \ U P 1 Vpl Upl' 

_ e— [5(r - r' - — )0(r' - r" - — ) + 0(r - r' - — )8{r' - r" - — )] .(A.* 

Upl Upl Upl Upl Upl 



Thus, the third-order contribution to Eq. I I A. 71) renomalizes the second-order 
term with 



M M 



L/v p i L/v p i 



2e E 9t>9t»1b9tCiM,\ f dr' f dr n 6(r' - r" - — ) x 

£'/"=l £ Upl 

a- 1+fc * , a- 1+fc ' , ^(T y )$ r (r") . (A.9) 
By setting e = ln(a/a ) in Eq. ( 1A.6I) . it is straightforward, but tedious, to derive 
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explicit expression for the nonlinear terms appearing in the RG equations for 
the running coupling strengths; the final result is 

,/;/ ' [< " [l-h e ]g e (a)+ E c £,£',e"9£'(a)9£"(a) . (A.10) 



d\n(a/ao) 



The RG equations in Eq. (1A.10l) may be rewritten as 
dgt dC[{g}] 



d\ia(a/a ) dgi 



where C[{g}] is given by 



, (A.ll) 



■y N y M 

c \{9}] = ^ El 1 - Hal + o E Ci,e',t"9e9i'9i" ■ ( A - 12 ) 
1 1=1 6 ii'i"=\ 



The RG fixed points coincide with the extrema of the function C[{g}], that 
is, with the set of values of {g}, {g*}, such that 

^M)l = , W=1,...,M . (A.13) 

dge 



Due to its properties, C may also be identified with the "boundary entropy" 
of the system [32] . 

We would like to point out that the thermodynamics limit may be achieved 
either by sending either a or L to oo. In fact, sending a to oo amounts to cut 
off high momenta contributions, which amounts to send the size of the system 
to oo. This implies that scaling may be realized either using a/a , or L/L* (as 
it has been done in this paper), as scaling parameters. 



B Boundary vertex operators 



As an example of correlation functions computed with both (3 and L finite, 
in this appendix we derive the correlators of the normal-ordered vertices 
Vx,±i(t) =: exp [±iX(ir)} :, and V^±i(ir) =: exp[±^(ir)] : introduced in 
section 3, where X(x,t) obeys Neumann boundary condions at both bound- 
ary x = and x = L, while t) obeys Dirichlet boundary conditions at the 
inner boundary, Neumann at the outer boundary. 



48 



The partition function Zx = Tr[— /3Hx], where Hx is the Hamiltonian for the 
field X, has been derived in Eg. (1841) . The result is 



Z X = 7~T eX P 



2gnv p i 2 
— p — — —m 



(B.l) 



The two- vertex correlation function, Gx{ti,t 2 ), is then given by 



1 

Z~x 



G x (r 1 ,r 2 ) = 7^Tr[T r {: e ~ iX(lri) :: e iX{iT ^ :}e~ fiHx ] 



(B.2) 



As it happens for the partition function, also Gx is factorized into a contri- 
bution from the oscillatory modes, times a contribution from the zero modes. 
The trace over the oscillatory modes yields 



IIT, 



- OSC 1 1 "*" n ' 

n=l 



(B-3) 



where 

T n = TrT r ^xp 

exp 



. Jla(-h) eVplknT1 



g n 



exp 



V g n 



j2ax(-k n ) eVplknT2 
V 9 n 



exp 



g n 



q 



ax{-kn)ax{kn) 



.(B.4) 



The trace in Eq. flB.4[) may be computed by resorting to a coherent state 
decomposition, which uses the basis of coherent states \x n ), defined by 



ax(k n )\x n ) = Vnx n \x n ) =>> \x n ) = e 2 e 



lire I 'n'xt-*") 



10} 



(B.5) 



and the decomposition of the identity given by 



(-) J dx n e ( 2' \x n )(x n \ 



(B.6) 



The derivation is tedious, but straightforward [37]. The final result is 
exp I I x 



1 



1 - q r > 



gn \ 1 — q r 
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cxp 



\ + / ■ \ 



M 1 - Q 



n 



n 



(B.7) 



As a consequence, one gets 

1 



IIT, 



A 00 

E 



X 



exp 



00 ^ e -fcnfpi|ri~r2|^?i -y 00 ^( e —knVpi\ri—Ta\\ri 

s£ (1 - 9") + </£ (1 - 9 n ) 



(B.8) 



Using the Jacobi's triple product identity |37J, one finally obtains 
00 1 rg a ^V(n-^) |^p_i^^ 



2L I L / 



n=l 



77(g) I 0; (0|i^) 



(B.9) 



where 9i((\t) is Jacobi's 6*i elliptic function, defined as 
O^Clr) = -iJ2 (-l) n e" r ( n+ ^) 2 e 2i7rC ( n+ ^) . 

neZ 



(B.10) 



By tracing over the zero mode and applying Poisson's resummation formula, 
one eventually gets 



Gx(n,T 2 ) = 9 3 



TU 



2gL 



. f3u s 



0'i (0|<^) 



(B.ll) 



where the ^-function is defined as 
3 (C\r)= E^^^ ■ 



(B.12) 



By performing similar computations, one obtains 
1 

Z x 



G 6 {t u t 2 ) = — — Tr [T T {: e~^ lTl) :: e i?(iT2) :}e _/JH e] 



TU 

2g~L 



'04 ( 



ZUpl 



2L 



01 (o|^) 



(B.13) 
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where the #4-function is defined as 



(B.14) 



When taking the limit (3 — > oo, the correlators in Eqs.( 1B.lTllB.13l) give the 
ones that we have used in the paper. 
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